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We present a method, the dynamical cumulant expansion, that allows to calculate quantum 
corrections for time-dependent quantities of interacting spin systems or single spins with anisotropy. 
This method is applied to the quantum-spin model H — —Hz{t)Sz + V{S) with _ffz(±co) — ±00 
and *(-oo) = |-5> we study the quantity P{t) = (l - {Sz)JS) /2. The case V{S) = ~H^S^ 
corresponds to the standard Landau-Zener-Stueckelberg model of tunneling at avoided-level crossing 
for N = 2S independent particles mapped onto a single-spin-S* problem, P{t) being the staying 
probability. Here the solution does not depend on 5* and follows, e.g., from the classical Landau- 
Lifshitz equation. A term —DS^ accounts for particles' interaction and it makes the model nonlinear 
and essentially quantum mechanical. The 1/S corrections obtained with our method are in a good 
accord with a full quantum-mechanical solution if the classical motion is regular, as for D > 0. 

PACS numbers: 03.65.-w, 75.10.Jm, 05.45.Mt 
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The classical limit of quantum spin systems has been 
discussed for decades in terms of spin coherent states 
Q IE Q 111' Yet little is known, in particular, about 
quantum corrections to the classical dynamics of systems 
with large spin S. One possible reason is that initially 
narrow quantum wave packets are spreading with time t, 
i.e., quantum states depend on the combined parameter 
of type / S {a > 1). Therefore the limits t 00 and 
5 — > 00 do not commute and a perturbation theory in 
pure 1/S is usually impossible. One is left here with the 
question of how well defined asymptotic {t — s- 00) quan- 
tities behave in the limit S 00. Are they continuous 
functions of l/S for 1/S ^0? 

One of the problems of this kind is a recent gener- 
alization Mid of the Landau-Zener-Stueckelberg (LZS) 
problem l^, 0| of transitions at avoided level cross- 
ings that is a well known quantum effect ubiquitous in 
physics of atomic and molecular collisions. Still the stan- 
dard LZS effect formally allows classical description since 
the Schrodinger equation for a pseudospin S ~ 1/2 is 
equivalent to the classical Landau-Lifshitz (LL) equation 



m = 7[mxHoff], B.cs = -dH/dm, 



(1) 



where m = (S) /S", 7 = g^^/h, and, in particular, Hcg = 
H = HxBx + Hz{t)ez with — A (level splitting) and 
H^{t) = E_i{t)-Ei\t) (level bias) satisfying Hz[±(x,) = 
±00. Typical time dependence of is linear, = vt, 
or nonlinear The probability to stay at the initially 
populated level —1 is 



Pit) = [1 - m^it)] /2, 
and for the linear sweep d, , 

P{oo) = P = e-', 



P{-oo) = 1 



2hv 
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This "degeneracy" between quantum and classical de- 
scriptions of the LZS effect disappears if tunneling parti- 
cles interact with each other, as is the case for molecular 
magnets ^3 (see Ref. ^| for a recent review) or for 



the Bose-Einstein condensate j3 0| . Finding the 
time-dependent wave function for the whole system of 
iV ^ 1 particles that is specified by 2^ coefficients is of 
course a tremendous problem. As a plausible first step 
one can consider a simplified model in which each parti- 
cle interacts with all TV — 1 other particles with the same 
strength J 0, that maps onto a large-spin model 
(we use gii^ = h = 1 below) 

H = -H,{t)Sz + V{S), S^N/2 (4) 
ViS) = -H,Sx-DSl D = 2J, (5) 

whereas P{t) is still defined by Eq. ((2Jl. On one hand, our 
original model of interacting particles |^, in the limit 
iV — > cx) is described exactly by the mean-field approx- 
imation (MFA) that simplifies the problem to a nonlin- 
ear Schrodinger equation for a single tunneling particle 
IE Q, |15,, 16] . On the other hand, the mean- field limit 
N ^ 00 corresponds to the classical limit S 00 oi 
Eq. Q that results into the LL equation Eq. with 
classical energy TC = —Hz{t)mz — w(m), where [not to be 
confused with the sweep rate v\ v{va) = —Hxnix — drnj, 
and d = SD. However for iV 7^ 00 that better suits to real 
systems with finite range interactions one is confronted 
with an essentially quantum effect. 

One can also look upon the large-spin problem of Eq. 
Q from another perspective and consider it as a kind of 
scattering problem for a spin: Sweeping the field across 
the region of strong interaction, Hz(±oo) — ±00, with 
the task to find (S'^)^^^^, if in the initial state the spin 
was down, {Sz)t^_^ac = '^^^^ the "Quantum 

Nonlinear Switching (QNS) model" . This model is more 
general than the LZS model and it can be applied to real 
large spins rather than to large pseudospins composed of 
many pseudospins 1/2. It could be relevant for molecular 
magnets, and one can even include in ^(S) terms that 
are absent in the LZS tunneling model, e.g., the biaxial 
anisotropy £^(5"^ — Sy) characteristic to Fes. 

The purpose of this Letter is to investigate the semi- 
classical limit S ^ 1 and to work out an analytical 
method to calculate 1/S corrections to the mean-field 
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result for the staying probability P for the QNS model. 
This provides a crucial test for the MFA that was ap- 
plied to the LZS effect in systems of interacting parti- 
cles in the absence of a more accurate and still tractable 
method^, 0- In comparison to statics, apphcation of 
the MFA to dynamics is questionable since the stabil- 
ity of the solution of the nonlinear LL equation Eq. (QJ 
with respect to quantum fluctuations is not guaranteed. 
It becomes clear as soon as one realizes that dynami- 
cal MFA describes classical dynamics that might exhibit 
chaos, that is, classically computed quantities P{t) might 
be poorly defined. 

The method we propose here is a dynamic cumulant 
expansion (DCE) based on the normalized and sym- 
metrized spin cumulants that are defined by 



TOq = {Sa)/S, 

1 ({S^Sp) + {SpSa,) 



map = ^ 



(Sa) {S„) (6) 



etc. with a,f3 — x,y,z = 1,2,3. We use the Heisen- 
berg representation for spin operators Sa and compute 
matrix elements with respect to the initial state that 
is the particular case 9 = n oi the spin coherent state 



Cm — 



25* 

S + m 



1/2 



with 



S+m 



COS ■ 



J{S-m)cp 



(7) 



P, 13, 0, 0- Matrix elements of products of operators 
factorize in the classical limit, and, correspondingly, nor- 
malized cumulants of n spin operators are as small as 
Thus a theory formulated in terms of cumu- 
lants yields a quasiclassical expansion. It is sufficient to 
retain pair cumulants and drop higher-order cumulants 
to obtain 1/5* corrections. If the state of the system 
would be close to a coherent state during all the time, 
the cumulant method would yield the result without any 
difficulties. We will see, however, that this assumption 
fails at large times because the nonequidistant spectrum 
of a nonlinear spin system leads to dephasing of differ- 
ent |?7i)-components of the coherent state. One can think 
of a strongly localized cloud of classical spins that, how- 
ever, are precessing with slightly different frequencies. 
Whatever narrow is this cloud, it will spread along the 
classical trajectory within a dephasing time Tdeph- That 
is, cumulants cease to be small at large times and the 
DCE diverges. Mathematically it means, as we will see 
below, that P{t) depends on both and on the com- 
bined parameter S~^t'^, where the pure describes the 
quantum correction to the QNS and S~^t^ describes the 
dephasing effect. The DCE picks both terms and thus it 
diverges at t — > c». Fortunately this divergence can be 
eliminated to yield a pure 1/5* correction to P = P{oo), 
as shown in Fig. ^for the model described by Eq. 

To construct DCE, we use the Heisenberg equation of 
motion for spin operators with an interaction V^(S) that 
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FIG. 1: Staying probability P{t) for the quantum nonlinear 
LZS effect, a - Exact numerical solutions for 5 = 5 and S —> 
oo; b - Cumulant solution of order 1/5": Raw and antiphased 
(see text). 



is taken to be a fully symmetrized function of Sa 

Sa — ^aPjH^^{t,S)Si3 

Hfit, S) = H^it) - dViS)/dSj. 



(8) 
(9) 



Eq. ||SJ) is in fact fully symmetrized but we don't write 
it explicitly to save space. Our quantity of interest is 
ma{t) — {-'S\Sa{t) \ — S) /S. Taking this matrix element 
makes Eq. nonclosed because dV{S)/dS.y is an op- 
erator. After expanding the matrix element of Eq. ^ 
up to the pair cumulants and dropping higher-order ones 
being of order 0(1/5^) it becomes an equation for 
coupled to the pair cumulants map- Equation for the lat- 
ter also follows from Eq. ||SJ| and it is closed at order 
1/S. Initial conditions to these equations at t = —oo are 
rriz = — 1 and irixx = fUyy = 1/ (25) while other compo- 
nents are zero. As equations for ijia and map are cou- 
pled, its solution generates all powers oil/ S that should 
be suppressed to obtain a pure 1/5* expansion. To this 
end, we expand 



(1) 



2S 



T-aP 
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It is convenient to express the results in terms of 

i;„,...„Jm) = 5"-V,,...„„(5m) (11) 



T4,...„„(S) = 



(12) 



As a result one obtains the three equations 

Ai") = T.,(i,M(''^)A^r (13) 
= T„,(t,/x("))A^i'j+T'^,(t,M("))A^W (14) 

/i« = f,^(t,/.(«))A.J,^)+fc„;3,(M(°))MS, (15) 

where 

np{t,m) = e„;37-ff'^(i,m) (16) 

^^^(t, m) = Taf3{t,ni) - €ajsm^vsp{ni) (17) 

with i?°^(t, m) ^ H-y{t) - v^(m) and 

ka^^im) = -eQ;55U5^(m) - eQ5^^i;;5-y^(m)TO5. (18) 

Equations p3|l -p5 () should be solved numerically in or- 
der of their appearance. The closed Eq. (|13() is a non- 
linear Landau-Lifshitz equation, whereas Eqs. H14() and 
()15|l are linear. Although Eqs. H13(l - (|15|l and the method 
of their derivation resemble decoupling schemes that are 
usually not based on small parameters, the DCE is a rig- 
orous (although formal) expansion in powers of 1/S'. 

Let us now check Eqs. (|13ll - (|15|) for the exactly solv- 
able toy model of Eqs. Q with = 0, starting at t = 
with an arbitrarily directed spin coherent state, Eq. l|7|l. 
Their solution is nothing else than a l/S* expansion of 
the exact solution taken for 



t < Tree = tt/L* = TrS'/d, 



(19) 



where Tree is the recurrence period, i.e., l/S* expansion 
of 



m:,(t)±imj,(t) =e=F^('^*-'^)exp ^ sin^ 6* sin6l. 

(20) 

In fact this expansion is only valid for t <^ Tdoph, where 



Tdoph = VS/d < Trt 



(21) 



is the dephasing time mentioned above Eq. 0. 

Consideration of the toy model shows that spin switch- 
ing and dephasing are different effects that are well sep- 
arated in time for 3> 1: QNS occurs in the vicinity of 
the resonace, Tqns S^, while dephasing occurs much 
later, Tdeph S^^"^. Thus one concludes that the diver- 
gence of DCE because of dephasing is harmless and it 
can be cured. In the QNS setup, dephasing affects dy- 
namics of mz(t) via T^(S). For S —^ oo, time dependence 
of niz far past the resonance can be found perturbatively 
in v{m) / Hz{t) <C 1. Here precession of m around the z 
axis is much faster than the temporal change of ff^, i.e.. 



Hz{t)/H^{t) <C 1. Thus the classical energy Ti is nearly 
conserved over the period of precession, and it yields 



(22) 
(23) 



where m'^''^(i) corresponds to the asymptotic form 
Jii-^^) = —Hz{t)mz. Precession of m^-^^\t) causes oscil- 
latory behavior of Smz{t) and thus of P{t) that however 
vanishes in the limit Hz{t) oo (see Fig. ^). For finite 
5^1, dephasing of Sm^it) follows from the first-order 
quantum-mechanical perturbation theory 



m,ra' — — S 



,{t) 



(24) 

where = {m\V{^)\m') , A$w(0 = *m(0 - 

$„'(<), and $„(t) = -m fdt'Hz{t') + V^n,t. Dephas- 
ing stems from the nonequidistant spectrum, i.e., from 
the nonlinear m-dependence of Vmm- As coefficients of 
the unperturbed state Cm after crossing the resonance 
are close to those of Eq. Q and are localized around a 
classical value toq — cos 9, calculation in Eq. (|24|l yields 
results similar to those for the toy model. Expanding 
Sniz {t) in 1/5 similarly to that of Eq. H20|l yields a lin- 
ear time divergence 



6mz{t) cc tS^^ cos [^H, (t) + ilj{t)] 



(25) 



for the linear sweep Hz{t) — vt, where ^H^{t) = 
fdt'Hzit') andi/'(t) is a slowly varying phase generated 
by the interaction: ■ip{t) ~ w(m(t)). This is exactly the 
raw result of DCE shown in Fig. 

We have found the source for the divergence of the 
DCE and seen that the divergent contribution dmz{t) in 
fact tends to zero, if described properly by Eq. (|24|) while 
the QNS effect with 1/S corrections is entirely contained 
in the asymptotic value 771^(00) = mi^^' in Eq. H22|l . 
Now it is clear that one should simply project out parts 
of mz{t) that are oscillating at large times. The latter 
can be done, e.g., by the "antiphasing" 



t + 



Thz 



t - 



Thz 



, (26) 



where Thz — '^'^/Hz{t) is the period of oscillations. This 
transformation kills the diverging factor t in Eq. (|25|) . To 
improve convergence, one can repeat antiphasing several 
times, including that with multiple precession frequencies 
arising due to nonlinearity. The result shown in Fig. 
at t ^ 00 is in a good accord with the exact numerical 
solution of the Schrodinger equation for P = P{oo) that 
is shown by a horizontal line. 

Fig. 121 shows P vs the quantum parameter 1/(25') for 
different sweep-rate parameters s of Eq. Q . The agree- 
ment between the exact numerical solution and the result 
from DCE for d = SD = 3 (that corresponds to strong 
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FIG. 2: Quantum corrections to the final staying probabifity 
P become more pronounced for slow sweep, e S> 1. 

ferromagnetic interaction) and = 1 is very good. Par- 
ticularly surprizing is the linearity of P vs 1/(25') down 
to S' = 5/2! Quantum corrections to P become more 
pronounced for slow sweep, £ ^ 1, and the slope of the 
curves tends to infinity in the limit e oo. This is in ac- 
cord with the finding of Ref. 7J: For = Hx/{2d) < 1 

in the classical case limg^oo -P = (1 ~ hl^^)'^^'^, whereas 
in the quantum case lime^oo P = because of high-order 
tunneling analogous to tunneling in molecular magnets. 

We have seen above that divergence of DCE because of 
dephasing is in fact a formal problem that can be solved 
to yield an excellent description of quantum effects in the 
QNS model. However there is another more fundamental 
problem related to the smoothness of the classical solu- 
tion. For the antiferromagnetic sign of the interaction, 
d < 0, the classical solution for P is a nonmonotonic 
function of e dropping fast to P = close to critical val- 
ues e,y, ly — 1,2, . . . (see Figs. 6 and 8 of Ref. 0)- 
are bifurcation points for classical trajectories at which 
for finite S the quantum distribution |cm(i)|^ that was 
a narrow packet at the beginning splits into two pack- 
ets when crossing the resonance. One can say that two 
parts of the quantum packet follow different classical tra- 
jectories. In this case DCE cannot be expected to work 
well since the assumption of a narrow packet close to a 
coherent state is not fulfilled. Accordingly quantum cor- 
rections are very large, and one needs an extremely large 
spin to approach the classical limit near these points (see 
Fig. 8 of Ref. 0). The curves P vs 1/(25) are not 
straight in the scale of Fig. |21and they are very different 
for different e. 

Classical models with more complicated interactions 
such as the biaxial model with the transverse field 
along the hard axis show dynamical chaos for time 
dependent Hz that makes the classical nonlinear spin 
switching effect poorly defined. In the chaotic regime the 
Schrodinger equation provides a more adequate descrip- 



tion of the problem than classical equations of motion 
as it "smoothes" the chaos (see, e.g., Ref. ^3). As a 
result, quantum effects are much stronger than 1/5' cor- 
rections in case of regular motion. On the other hand, 
for the field along the medium axis classical motion is 
regular and smoothly depends on parameters, thus DCE 
provides an excellent description of quantum corrections 
at order 1/5' that is comparable with that of Fig. El 

In this Letter, we have tested the dynamical cumulant 
expansion for the particular model of interacting tunnel- 
ing species with a special interaction that is equivalent to 
a large spin. We have shown that in the case of a smooth 
classical motion it works well. This method can be gener- 
alized for systems with realistic interactions that cannot 
be mapped onto a large spin and thus are problematic if 
one uses the Schrodinger equation for the whole system. 
The system of cumulant equations will include that for 
the magnetization nia^i that in the simplest case is in- 
dependent of the site i and those for the pair cumulants 
rriap^ij that depend on the distance between i and j. One 
has to check however before applying DCE that the clas- 
sical (i.e., the mean- field) solution of the problem is not 
chaotic. 
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